Measurement of K^ 3 form factors 



NA48 Collaboration 
A. Lai, D. Marras 

' Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Cagliari, 

O ■ 1-09100 Cagliari, Italy 

(N 

A. Bevan, R.S. Dosanjh 1 , T.J. Gershon 2 , B. Hay, G.E. Kalmus, 
C. Lazzeroni, D.J. Munday, E. Olaiya 3 , M.A. Parker, 
T.O. White, S.A. Wotton 



Cavendish Laboratory, University of Cambridge, Cambridge, CBS ORE, U.K.' 

G. Barr 5 , G. Bocquet, A. Ceccucci, T. Cuhadar-Donszelmann 6 , 
D. Cundy 7 , G. D'Agostini, N. Doble 8 , V. Falaleev, 
r^; ; L. Gatignon, A. Gonidec, B. Gorini, G. Govi, P. Grafstrom, 

W. Kubischta, A. Lacourt, A. Norton, S. Palestini, 



O 



V : B. Panzer-Steindel, H. Taureg, M. Velasco 9 , H. Wahl 1 ' 

43 



D . CERN, CH-1211 Geneve 23, Switzerland 



X 



C. Cheshkov, u ' 12 , P. Hristov, n ' 12 , V. Kekelidze, L. Litov 12 , 
D. Madigozhin, N. Molokanova, Yu. Potrebenikov, 
S. Stoynev, 9 ' 12 , A. Zinchenko 

Joint Institute for Nuclear Research, Dubna, 141980, Russian Federation 

I. Knowles, V. Martin 9 , R. Sacco 13 , A. Walker 

Department of Physics and Astronomy, University of Edinburgh, JCMB King's 
Buildings, May field Road, Edinburgh, EH9 3JZ, U.K. 



Preprint submitted to Physics Letters B 



3 February 2008 



M. Contalbrigo, P. Dalpiaz, J. Duclos, P.L. Frabetti 14 , 
A. Gianoli, M. Martini, F. Petrucci, M. Savrie 

Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Ferrara, 1-44100 

Ferrara, Italy 

A. Bizzeti 15 , M. Calvetti, G. Collazuol 8 , G. Graziani, 
E. Iacopini, M. Lenti, G. Ruggiero, M. Veltri * 16 

Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Firenze, 

1-50125 Firenze, Italy 

H.G. Becker, K. Eppard, M. Eppard 11 , H. Fox 17 , A. Kalter, 
K. Kleinknecht, U. Koch, L. Kopke, P. Lopes da Silva, 

P. Marouelli, I. Pellmann 18 , A. Peters 11 , B. Renk, 
S.A. Schmidt, V. Schonharting, Y. Schue, R. Wanke, 
A. Winhart, M. Wittgen 19 

Institut fur Physik, Universitdt Mainz, D-55099 Mainz, Germany 20 

J.C. Chollet, L. Fayard, L. Iconomidou-Fayard, J. Ocariz, 
G. Unal, I. Wingerter-Seez 

Laboratoire de I'Accelerateur Lineaire, IN2P3-CNRS,Universite de Paris-Sud, 

91898 Or say, France 21 

G. Anzivino, P. Cenci, E. Imbergamo, P. Lubrano, 
A. Mestvirishvili, A. Nappi, M. Pepe, M. Piccini 

Dipartimento di Fisica dell'Universita e Sezione dell'INFN di Perugia, 

1-06100 Perugia, Italy 

L. Bertanza, R. Carosi, R. Casali, C. Cerri, M. Chilli 11 , 
F. Costantini, R. Fantechi, S. Giudici, I. Mannelli, 
G. Pierazzini, M. Sozzi 

Dipartimento di Fisica, Scuola Normale Superiore e Sezione dell'INFN di Pisa, 

1-56100 Pisa, Italy 

J.B. Cheze, J. Cogan, M. De Beer, P. Debu, A. Formica, 
R. Granier de Cassagnac 22 , E. Mazzucato, B. Peyaud, 
R. Turlay, B. Vallage 

DSM/DAPNIA - CPA Saclay, F-91191 Gif-sur-Yvette, France 



2 



R. Bernhard 23 , M. Holder, A. Maier 11 , M. Ziolkowski 

Fachbereich Physik, Universitdt Siegen, D-57068 Siegen, Germany 24 " 

R. Arcidiacono, C. Biino, N. Cartiglia, F. Marchetto, 
E. Menichetti, N. Pastrone 

Dipartimento di Fisica Sperimentale deU'Universitd e Sezione dell'INFN di 

Torino, 1-10125 Torino, Italy 

J. Nassalski, E. Rondio, M. Szleper 9 , W. Wislicki, S. Wronka 

Soltan Institute for Nuclear Studies, Laboratory for High Energy Physics, 
PL-00-681 Warsaw, Poland 25 

H. Dibon, G. Fischer, M. Jeitler, M. Markytan, I. Mikulec, 
G. Neuhofer, M. Pernicka, A. Taurok, L. Widhalm 

Osterreichische Akademie der Wissenschaften, Institut fur Hochenergiephysik, 

A- 1050 Wien, Austria 26 



Abstract 

This paper reports on a new high precision measurement of the form factors of the 
Kl — ► tt^ /j/^ Va decay. The data sample of about 2.3xl0 6 events was recorded 
in 1999 by the NA48 experiment at CERN. Studying the Dalitz plot density we 
measured a linear, A + = (20.5±2.2 st at ±2.4 S2/st ) x 10 -3 , and a quadratic, A + = (2.6± 
0.9 s tat i 1.0 S yst) x 10 -3 term in the power expansion of the vector form factor. No 
evidence was found for a second order term for the scalar form factor; the linear slope 
was determined to be Ao = (9.5±l.l siai ±0.8 sys t) x 10 -3 . Using a linear fit our results 
were: A+ = (26.7±0.6 stat ±0.8 syst ) x 10" 3 and, A = (11.7+0.7,^1.0^) x lO" 3 . 
A pole fit of the form factors yields: my = (905 ± 9 s t a t ± 17,^) MeV/c 2 and 
m s = (1400 ± 46 stat ± 53 syst ) MeV/c 2 . 



Key words: Kaon semileptonic decays, Kaon form factors 
PACS: 13.20 Eb 
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1 Introduction 



Since long ago [1] K £3 decays (K L — > ir ± £ T v, £ = e, /x) have offered the 
opportunity to test several features of the electroweak interactions such as 
the V-A structure of weak currents, current algebra and the predictions of 
chiral perturbation theory. These decays have been the object of a renewed 
interest both on the experimental and theoretical side since they provide the 
cleanest [2] way to extract the CKM matrix element \V US \. K £3 decays give 
access to the product / + (0)|V^|, where /+(0), the vector form factor at zero 
momentum transfer, has to be determined by theory. The recent calculations at 
0(p 6 ) [3] in the framework of chiral perturbation theory show how /+(0) could 
be experimentally constrained from the slope and the curvature of the scalar 
form factor / of the K^ 3 decay. In addition, the form factors are needed 
to calculate the phase space integrals which are another ingredient for the 
determination of \V US \. Finally, very recently it has been pointed out [4] how 
a precise measurement of the value of fo at the Callan-Treiman point [5] 
could provide a clean test of a small admixture of right handed quark currents 
(RHCs) coupled to the standard W boson. 

Until recently, the experimental knowledge [6] on K^ 3 form factors was mainly 
based on a certain number of old measurements dating back to the seventies. 
The slopes obtained from K M 3 decays were less precise than those determined 
in K C 3 decays, and a large difference between the results from charged and 
neutral kaon decays was present. This difference was more pronounced for the 
slope Ao where, in addition, the situation was confused also by the presence of 
negative values. Very recent high precision experiments [7,8,9,10,11] provided 
a more accurate determination of these quantities with values smaller than 
the old PDG averages and agreement between K° and measurements has 
been established. Furthermore, evidence for a quadratic term in the vector 
form factor was found, at the level of about 2 a, by ISTRA+ in and by 
KLOE in K e3 decays. A cleaner indication came also from KTeV, both in 
K e3 and decays, with a significance of about 3 a. 

This paper reports on a new high statistics measurement of K^3 form factors. 
Following this introduction Sec. 2 recalls the formalism about the decays, 
Sec. 3 describes the experimental set-up, Sec. 4 reports about the analysis, 
and Sec. 5 delineates the fitting procedure and the treatment of the systematic 
error. 



using computing resources of the Interdisciplinary Center for Mathematical and 
Computational Modelling of the University of Warsaw. 

26 Funded by the Federal Ministry of Science and Transportation under contract 
GZ 616.360/2-IV GZ 616. 363/2- VIII, and by the Austrian Science Foundation under 
contract P08929-PHY. 
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2 The K L -> 7r ^ ^ decay 



Only the vector current contributes to K^ 3 decays. As a result the matrix 
element can be written in terms of two dimensionless form factors f± : 

Tl=-^=V US [/+(*) (P K + P^f uty^l + 75K + f-(t) m t u,{\ + lh )u v \ , (1) 

where Pk and are the kaon and pion four momenta, respectively, ui and 
u v are the lepton spinors, rri£ is the lepton mass and t = (P# — P-*) 2 = 
m 2 K + m 2 — 2Pk ■ P n = q 2 is the square of the four-momentum transfer to the 
lepton pair. The form factor /_ (t) is related to a scalar term proportional to 
the lepton mass and can be measured only in K m3 decays. 
The determination of the form factors in this analysis is based on a study of 
the Dalitz plot density which can be parametrized [1] as: 

p(e;, k) = dN d ^E* :) K Af + {t) + B f+W-® + Cf2 ^ ( 2 ) 

where A, B and C are kinematical terms: 

A = 7jik(2E*E* - m K El) + (1/4 < - 
B = ml(E;- 1/2 E' n ), 
C = 1/4 mlK- 

E*, E* are the muon and pion energy in the kaon center of mass (CMS) re- 
spectively. For the neutrino we have E* = rrtx — E* — E* and E' w is defined 
as: 

Tp rp*Max rp* _ n M rp* 

n ~ * 77 ~ 2m K 

In an alternative parametrization one can define the form factor /o(^) : 

/oft) = /+(*) + 7 2 - 2 , /-(*)• (3) 

i m K ~ m J 

This parametrization is preferred because /+ and f are related to the vector 
(1~) and scalar (0 + ) exchange to the lepton system, respectively and are less 
correlated than in the previous case. The expansion in powers of t of the form 
factors is often stopped at the linear term: 

/±,o(*) = /±,o(0) (l + \ ± , t/miy (4) 

The assumption that /+ and / are linear in t implies that /+(0) = /o(0) so 
that /_ does not diverge at t = 0. 

The form (4) is usually adopted as consequence of the smallness of data sam- 
ples from the past experiments rather than on physical motivations. Nowadays, 
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with higher statistics, it is becoming customary to search for a second order 
term in the form factors expansion: 



/+,„(*) = /+(0) 



1 + A o t/< + ^ Al {t/mlf 



(5) 



The weak interaction of hadron systems at low energies can also be described 
in terms of couplings of mesons to the weak gauge bosons (pole model, meson 
dominance [12]). In this framework the form factors acquire a physical mean- 
ing since they can be related to the exchange of the lightest K* resonances 
which have spin-parity l~/0 + and mass m v /m s , respectively: 



/+(*) = /+(0) 



m 



v 



t 



fo(t) = /+(0) 



m\ 



t 



(6) 

o v - l m" s -i 

Recently new parametrizations of the vector [13] and scalar [4] form factors 
based on dispersion relations subtracted twice have been proposed: 



/+(*) = /+(0)exp 



t 



— (A ++ H(t)) 



rat 



/ o (i) = / + (0)exp 



t 



[m 



A 



mi 



■(lnC-G(t)) 



(7) 

here A + is a slope parameter and InC = \n.[fo(m 2 K — m*)] is the logarithm 
of the value of the scalar form factor at the Callan-Treiman point. For the 
dispersive integrals H(t) and G(t) accurate polynomial approximations have 
been derived. 



3 Experimental set— up 



For the measurement reported here the data were taken during a dedicated 
run period in September 1999. A pure Kl beam was produced by 450 GeV/c 
protons from the CERN SPS hitting a beryllium target. The decay region was 
contained in a 90 m long evacuated tube and was located 126 m downstream 
the target after the last of three collimators. 

The NA48 detector was originally designed for a precise measurement of direct 
CP violation in the neutral kaon decays to two pions. A detailed description 
can be found in [14]; only the main components relevant for this measurement 
are described here: 

Magnetic Spectrometer. It was contained in a helium tank kept at at- 
mospheric pressure and consisted of four drift chambers and a magnet. Each 
chamber had four views (x, y, u, v) each of which had two planes of sense wires. 
The spatial resolution per projection was 100 /im and the time resolution was 
0.7 ns. The magnet, placed between the second and the third chamber, was 
a dipole with a transverse momentum kick of 265 MeV/c. The momentum 
resolution of the spectrometer was (p in GeV/c): 
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= 0.48% © 0.009 p%. 

V 

Hodoscope. Located downstream of the spectrometer, it was used to pro- 
vide a precise time reference for tracks. It consisted of two orthogonal planes 
of scintillators segmented in horizontal and vertical strips and arranged in four 
quadrants. The time resolution per track was about 200 ps. The coincidence 
of signals from quadrants was used in the first level trigger for events with 
charged particles. 

Electromagnetic calorimeter. This was a quasi-homogeneous liquid kryp- 
ton calorimeter (LKr) with projective tower structure made by Be-Cu 40 jum 
thick ribbons extending from the front to the back of the device in a small 
accordion geometry. The 13248 read-out cells had a cross-section of 2 x 2 cm 2 . 
The energy resolution was parametrized as (E in GeV): 

a E (3.2±0.2)%„(9 ± l)% ffi(a42±a05)% 



E y/E E 

Muon system. The muon system (MUV) was located between the hadron 
calorimeter and the beam dump and consisted of three planes of scintillators 
each shielded by a 80 cm thick iron wall. The first two planes were made of 
25 cm wide horizontal and vertical scintillators strips. The strips overlapped 
slightly in order to ensure full coverage over the whole area of 2.7 x 2.7 m 2 . 
The third plane had horizontal strips 44.6 cm wide. The central strip was 
split with a gap of 21 cm to accomodate the beam pipe. All counters (apart 
the central ones) were read out at both sides. The inefficiency of the system 
was at the level of one per mill and the time resolution was below 1 ns. The 
passage of particles in the MUV produces "hits", i.e. a coincidence between 
an horizontal and a vertical counter which defines a 25 x 25 cm 2 region. 
Trigger. The acquisition of events was driven by a two level trigger. In the 
first level the presence of at least two hits in the hodoscope was requested. 
In the second level trigger the spectrometer was used to reconstruct tracks 
and a vertex made of opposite charge tracks in the decay region was required. 
To measure the trigger efficiency, a control trigger was implemented using the 
first level trigger properly downscaled. 



4 Data analysis 



4-1 Event selection 



The data sample consists of about 10 8 triggers recorded alternating the po- 
larities of the magnetic field of the spectrometer. To identify the K m3 decays 
the following selection criteria were applied to the reconstructed data. 
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The events were required to have exactly two tracks of opposite charge form- 
ing a vertex in the decay region, defined to be between 7.5 m and 33.5 m from 
the exit of the final collimator and within 2.5 cm from the beam line. The 
distance of closest approach of these two tracks had to be less than 2 cm. 
The difference in the track times reconstructed by the spectrometer had to be 
less than 6 ns while for the times determined by the hodoscope a maximum 
difference of 2 ns was admitted. 

Both tracks were required to be inside the detector acceptance by demanding 
that their projection had to be inside the fiducial area of the various subde- 
tectors. Tracks were accepted in the momentum range between 10 and 170 
GeV/c. 

In order to allow a clear separation of showers, a minimum distance of 35 
cm between the extrapolated impact points of the tracks at the entrance of 
the LKr calorimeter was required. Furthermore to avoid problems due to the 
misreconstruction of the shower energy a minimum distance of 2 cm from the 
track impact point to a dead calorimeter cell was imposed. 
In order to reduce the background from Kl — > n + n~ ir° (K 37r ) decays the 
cut P' Q 2 < -0.004 (GeV/c) 2 was applied. The variable Pq 2 , which is computed 
assuming that the decay is a K 3?r , is defined as: 

I 2 2 2 \ 2 A I 2 2 , 2 2 \ 

[m K - m + _ - m n o J - 4 [m + _m 7r o + m K p ± J . 

In the above formula p± (m + _) is the transverse momentum (invariant mass) 
of the assumed tt + tc~ system relative to the direction of the Kl- P 2 represents 
the Kl momentum in a reference frame in which the longitudinal component 
of the pion system is zero. It is positive for K 37r decays and negative for rQ 3 de- 
cays. 

The K e3 background was suppressed using the ratio E/p, where E is the energy 
of the cluster, reconstructed in the electromagnetic calorimeter and associated 
to a track, and p is the track momentum as measured by the spectrometer. 
For both tracks E/p had to be less than 0.9. The probability for a it to be 
rejected by this cut was measured to be about 1%. 

A track was identified as a muon when its extrapolated impact point at the 
MUV could be spatially associated to a MUV hit. The distance of association 
was dependent on the momentum of the track to account for multiple scatter- 
ing and measurement errors. For this analysis in addition, other constraints 
were imposed: the distance between the track extrapolation and the hit had 
to be less than 30 cm; the difference between the event time (determined by 
the charged hodoscope) and the muon time (determined by the MUV) had 
to be less than 3 ns, and finally also a coincidence in the MUV plane 3 was 
required. Monte-Carlo (MC) studies indicate that under these conditions the 
probability to misidentify a /i for a ir is at the 10~ 5 level. 
A well-known problem with PQ 3 decays is the quadratic ambiguity in the de- 
termination of the Kl energy. The v escapes undetected in this decay and 
while the transverse component of the momentum {p u t) to the Kl direction 



p'2 



4(v 2 



m 
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of flight (obtained joining the event vertex to the target position) is deter- 
mined by the /x and it transverse momenta, the longitudinal component (pIl) 
can be determined only up to a sign representing the two possible orientations 
of the v in the kaon CMS. This ambiguity leads to two solutions for kaon 
energy, called "low" (El) and "high" (Eh)- As an additional selection criteria 
we required both kaon energy solutions to be greater than 70 GeV. 
Finally a cut was applied on the variable p* — p U T, where is the total neu- 
trino momentum in the kaon CMS. This quantity, clearly positive for good 
K M 3 events, is highly sensitive to resolution effects which give rise to a mod- 
erate negative tail. We set a cut at p* — p v r > 7 MeV/c, selecting a region 
where the MC simulation accurately reproduces the data behaviour. After the 
applications of all cuts, 2344382 K m3 events were reconstructed from the data 
sample. 

x 10 2 



5000 
!T> 4500 



1 4000 



3500 
3000 
2500 
2000 
1500 
1000 
500 




i.i 

1.075 
1.05 
1.025 ~r 
1 

0.975 7 
0.95 ~r 
0.925 ]■ 
0.9 



Data/MC 



X/dof = 31.6/28 



80 100 120 140 160 180 200 220 



+ Data 

o MC 



75 100 125 150 175 200 225 250 275 



E, (GeV) 



Fig. 1. Comparison data-MC for the low kaon energy solution, the inset shows the 
ratio data/MC. 



4-2 Monte-Carlo simulation 



The detector response has been simulated in details using a MC program 
based on GEANT [15]. Particle interactions in the detector material as well as 
response functions of the different detector elements have been taken into ac- 
count in the simulation. Pre-generated shower libraries for photons, electrons 
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and charged pions are used to describe the response of the calorimeters. 
To determine the detector acceptance as well as the distortion and losses of 
events on the Dalitz plot induced by the radiative effects, the K m3 decay has 
been simulated both at the Born level and at the next-to-leading order (NLO). 
The acceptance suffers only from a residual dependence on the values (and on 
the type of parametrization) of the form factors used for the generation of 
the MC samples. To avoid any biases the samples were produced, after an 
iterative procedure, with values close to the results reported here. A linear 
parametrization was used with A+ n = 0.0260 and Aq 6 " - = 0.0120. 
The NLO sample was obtained using as event generator KL0R [16], a pro- 
gram which numerically evaluates the radiative corrections and generates MC 
events. The simulated events underwent the same reconstruction procedure as 
the data events and the same selection cuts described in Sec 4.1 were applied. 
These two MC samples provide a statistics which is one order of magnitude 
larger than the data one. 

A third, smaller, K^ 3 sample was generated with full simulation of the show- 
ers in the calorimeters and was used to model the multiple scattering in 
the MUV. For detailed studies of the background samples of K e3 , K 37r and 
K L — ► 7r + 7r~ (K 27r ) decays decays were produced. 

The Kl energy spectrum was extracted from the data by using the data dis- 
tributions of the low and high energy solutions and the probability matrix, 
obtained from MC, which relates El and Eh to the true kaon energy. 
To show the quality of the MC simulation the comparison of data and MC for 
some relevant kinematical quantities are shown in Fig. 1 to 3. 

4-3 Backgrounds 

The K e3 , K 3n and K 27r decays are the major sources of background. A K e3 event 
can fake a genuine K m3 when the 7r decays into a muon and the e has the E/p 
requested for the identification of a n (E/p < 0.9). To determine this source 
of contamination in the selected K m3 sample, we used the E/p distributions 
of tracks which pass all cuts, but not considering E/p. The electron signal is 
obtained by fitting this distribution around the value of 1. The integration of 
the fitted function into the "pion" region allows to determine a value for the 
K e3 contamination of: 

-p%£ = (6.59 ± 0.09) x 10" 4 . (8) 

The K 37r decays (followed by the decay of one of the two charged it) are 
strongly suppressed by the P Q 2 cut. To determine the residual contamination 
the selected K m3 events undergo a K 3n selection procedure: in the presence of 
clusters in the LKr not associated to the tracks, an attempt to reconstruct a 
7T° is made. In case the two photons reconstruct the it mass within a window 
of ±6 MeV/c 2 , the invariant mass of the two tracks (assumed to be pions) 
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Fig. 2. Comparison data-MC for the high kaon energy solution, the inset shows the 
ratio data/MC. 

and the tt° is evaluated and if it falls in an interval of ±9 MeV/c 2 around 
the K° mass the event is assumed to be a K 37r decay. The number of these 
K 37r background events, corrected for their acceptance, allows to estimate for 
the K 3?r contamination the value: 



1~)Cont 
1 3vr 



(6.31 ±0.16) x 10" 



(9) 



Another source of background stems from the decay with subsequent 7r 
decay in flight or pion punch-through in the iron of the MUV. Using the 
K 2n MC sample this contamination is estimated to be: 



'rycont 
1 2tt 



= (5.63 ±0.16) x 10" 4 . 



(10) 



This background source turns out to be the most dangerous one since the 
K 2n events populate a narrow region (the top right corner) of the Dalitz 
plot introducing appreciable distortions. The K 37r events instead populate the 
bottom left region of the plot; being not much concentrated, they induce a 
smaller effect. Finally the K C 3 events are distributed randomly on the plot and 
their effect is negligible. 

Background events from K 27r and K 37r will be subtracted from the data while 
the effects of K c3 events will be included in the treatment of the systematic 
uncertainty related to the background (Sec. 5. 2). 



12 




20 40 60 80 100 120 




20 40 60 80 100 120 



P % (GeV/c) 




0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 
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Fig. 3. Comparison data-MC for the momentum distributions of muons, pions (lab 
system) and neutrinos (kaon CMS). 

5 Fitting procedure and results 

5. 1 Fitting procedure 

The measurement reported here is based on the study of the Dalitz plot den- 
sity. As mentioned before, the ambiguity in the determination of the kaon 
energy leads to two solutions for the Kl energy and for the CMS energies of 
the /i and the n. Consequently each event has a double location on the Dalitz 
plot. We chose to evaluate E* and E* by using only the low kaon energy solu- 
tion. According to the MC simulation, this corresponds to the most probable 
solution, being in 61% of cases the correct one. The Dalitz plot was divided 



13 



I 0.26 



0.24 
0.22 
0.2 
0.18 
0.16 
0.14 









1 






















1 






! 








- 










































- 


















'1"i'"'i'"'s> 


















- 


































K 
































































































- 






































; 


— 


- 




































K 




"- 


































































TV 


i 




























w 






















— 






































— 


- 










































"■ 












































































































































- 
















































- 


- 














































"- 


- 
















































"- 
















. 


... 












































£> 
















































































































































































































8 




























































































































































































































































































































































































r 














tt 









80000 
70000 
60000 
50000 
40000 
30000 
20000 
10000 




0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 

E*\l (GeV) 

Fig. 4. Dalitz plot distribution of the K^3 events corrected for radiative effects and 
acceptance; the shaded cells are not used for the fit. 



into cells with a dimension of about 4x4 MeV 2 (see Fig. 4). About 39% of the 
events are reconstructed exactly in the same cell where they were generated, 
while this figure drops to 22% if the high solution is used. To extract the form 
factors we fit the data Dalitz plot, corrected for acceptance and radiative ef- 
fects, to the Born level prediction. The acceptance, in the z-th cell of the plot, 
6i, is defined as the ratio of the number of reconstructed events (evaluated 
using the low energy solution) to the number of generated events (evaluated 
using the true kaon energy) in that cell. We note that this definition of accep- 
tance accounts also for the migration of events induced by the use of the low 
solution only. 

The correction (for the i-th cell of the plot) due to the radiative effects is 
(1 + 5)i and is evaluated by taking the ratio between the number of recon- 
structed events from the MC-NLO sample and the number of reconstructed 
events from the MC-Born one. 

The number of events, corrected for acceptance and radiative effects, in a given 
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cell of the plot is therefore: 



Ni= * f , , (11) 



where A^ i Rec is the number of reconstructed and background subtracted data 
events. 

The form factors were determined by fitting with the MINUIT [17] package the 
Dalitz plot distribution, corrected for acceptance and radiative effects, to the 
parametrization reported in Eq. 2. The cells crossed by the Dalitz plot bound- 
ary are excluded from the fit (see Fig. 4). Various t dependences of the form 
factors were considered: linear, quadratic, pole and dispersive. The fit results 
are listed in Table 1; the correlations among the fitted form factors parameters 
are shown in Table 2. The comparison Data-Fit are shown in Fig. 5. 
We also fitted for a possible quadratic term in the scalar form factor and 



Linear (xlO 3 ) 


A+ 


Ao 


X 2 / n df 




26.7±0.6 


11.7±0.7 


604.0/582 


Quadratic (xl0~ 3 ) 


a; 


a; 


Ao X 2 /ndf 




20.5±2.2 


2.6±0.9 


9.5±1.1 595.9/581 


Pole (MeV/c 2 ) 


my 


m s 


X 2 /ridf 




905±9 


1400±46 


596.7/582 


Dispersive (xlO -3 ) 


A+ 


lnC 


X 2 /ndf 




23.3±0.5 


143.8±8.0 


595.0/582 



Table 1 

Form factors fit results for linear, quadratic pole and dispersive parametrizations. 
The quoted errors are the statistical ones. 

found Aq = (1.1 ± 1.3) x 10 -3 indicating that the linear assumption is suffi- 
cient to describe this form factor. In the expansion of the vector form factor 
instead, evidence is present for the existence of both a linear and a quadratic 
term. We notice also a remarkable shift in the value of Ao as consequence of 
the presence of the quadratic term in the vector form factor expansion. The 
value of my and ms, obtained with the pole fit are found to be consistent 
with the K*(892) and K*(1430) masses, respectively. 

As a cross-check we extracted the linear form factors using a x 2 built by com- 
paring the data Dalitz plot distribution, corrected for radiative effects only, 
with a set of Born level plots of reconstructed MC events. Each MC Dalitz 
plot distribution was produced with different form factors values by proper 
re-weighting of the events of the reference MC-Born sample. The form fac- 
tors values are extracted by minimizing the x 2 (A + , Ao) function constructed 
in this way. The results obtained with this method are less accurate than those 
provided by the fit procedure but fully unbiased with respect to the choice of 
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Linear 


A+ 


Ao 

n An 
-U.4U 




Quadratic 




a; 


Ao 




A + 


-u.yo 


U.Oo 




a; 




-0.73 


r ore 


my 


ms 
-0.47 




Dispersive 


InC 


A+ 
-0.44 





Table 2 

Correlation coefficients among the fit parameters for the linear, quadratic, pole and 
dispersive parametrizations. 

the form factors values used to generate the MC sample. These results are in 
perfect agreement with the ones obtained with the fit procedure, indicating 
the absence of such kind of bias in the analysis procedure. 
To check the fit procedure we fitted MC events, using the reference MC sample 
(generated with linear parametrization) and smaller samples generated with 
quadratic and pole parametrizations. In all the three cases the input form 
factors were correctly reproduced at the end of the process, indicating the 
absence of any bias in the fit procedure. 



5.2 Systematic uncertainties 

Various sources of systematic uncertainties in the determination of the form 
factors have been investigated. Their individual contributions are reported on 
Table 3 together with the effects related to the background contamination. 
The total error was obtained by combining the individual errors in quadra- 
ture. 

Effects related to the background have been checked altering the estimated 
contaminations by 15% and accounting for the tiny effect related to K e 3 events. 
The variations in the fit results were taken as the systematic uncertainty. 
Effects related to the acceptance and selection criteria have been checked by 
varying the selection cuts in a reasonable range. The largest fluctuations in 
the form factors were taken as systematic errors. 

Effects related to the Kl energy spectrum used in the MC simulations were 
investigated by using the spectrum obtained from a clean sample of K 27r de- 
cays. The simulated events were re-weighted with the ratio of the two spectra, 
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Fig. 5. Comparison Data-Fits and Data/Fit ratios for the four parametrization used 
in the analysis. For visualization purposes the various Dalitz plots distributions have 
been projected onto the E* axis. 



and the differences in the form factor results were taken as the systematic un- 
certainty. 

To check effects related to the use of the low kaon energy solution the analysis 
was repeated using the high solution to determine the acceptance and radia- 
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AA+ AA AA' + AA+ AA Amy Am s AA+ AlnC 
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1 
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±0, 
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Total Error ±1.0 ±1.2 ±3.3 ±1.3 ±1.4 ±19 ±70 ±0.9 ±13.8 

Table 3 

Systematic and total uncertainties for the four form factor parametrizations ana- 
lyzed. The systematic and statistical uncertainties have been added in quadrature to 
obtain the total error. 

tive corrections. Also in this case the differences with the reference fit results 
were taken as systematic error. 

The inefficiency of the MUV during this run was measured by identifying 
the /i according to its energy deposition in the electromagnetic and hadronic 
(HAC) calorimeters. The MUV efficiency was found to vary between 0.97 
for a 10 GeV/c and 1 for a > 20 GeV/c muon with an average of €muv = 
0.9987 ± 0.0001. To investigate possible biases, the inefficiency was artificially 
increased by randomly rejecting events according to the momentum depen- 
dence of the efficiency and its value, without observing any significant effect. 
Effects related to the MUV offline reconstruction were tested by relaxing the 
cut between the track extrapolation and the hit position and by accepting also 
events for which only plane 1 and 2, but not plane 3, had fired. This produces 
an increase of 3.2% in the statistics of the data sample; here again differences 
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from the reference form factor values were taken as systematic error. 
Effects related to the radiative corrections model used in the analysis were 
tested by applying the corrections obtained with the Ginsberg [18] formalism, 
amending the error reported in Ref. [19] and allowing for a t dependence of 
the form factors. The differences with the reference results were taken as an 
estimate of the systematic effect. 

The effects related to the size of the cells in which the Dalitz plot was divided 
were determind by reducing the cell size down to about 3x3 MeV 2 , the largest 
fluctuations in the form factors were taken as systematic errors. 
To estimate the possible influence of accidental particles, tracks outside the 
allowed time window for a match in the MUV were studied. No effect was 
found from this source. 



6 Conclusions 

The K^ 3 decay has been studied with the NA48 detector. A sample of 2.3xl0 6 
reconstructed events was analyzed in order to extract the decay form factors. 
Studying the Dalitz plot density we measured the following values for the form 
factors parameters: A' + = (20.5 ± 2.2 stat ± 2A syst ) x 10~ 3 , A" = (2.6 ± 0.9 sta t ± 
1.0 syst ) x 1(T 3 and A = (9.5 ± l.l stat ± 0.8 syst ) x 10~ 3 . 

Our results indicate the presence of a quadratic term in the expansion of the 
vector form factor in agreement with other recent analyses of kaon semileptonic 
decays. Fig. 6 shows the comparison between the results of the quadratic fits 
as reported by the recent experiments [11,9,10,7,8]. The 1 a contour plots are 
shown, both for K e 3 and K M 3 decays; the ISTRA+ results have been multiplied 
by the ratio (m n + /m^o) 2 . The results are higly correlated, those from this 
measurement and from KTeV have a larger quadratic term and appear only in 
partial agreement with the other K e3 experiments. We notice however that the 
observed spread in the A + , A + figures is greatly reduced if the values obtained 
from the Taylor expansion of the pole parametrization (A + = m^/my; A + = 
2 A. 2 ) are used. 

Using a linear fit our results were: A + = (26.7 ± 0.6 stet ± 0.8 syst ) x 10~ 3 and, 
A = (11.7 ±0.7 stat ±1.0 sy st ) x 10~ 3 . While the result for A + is well compatible 
with the recent (and most precise) KTeV measurement, the value of Ao appears 
to be shifted towards lower values. A pole fit of the form factors yields: my = 
(905 ± 9 stat ± 17 syst ) MeV/c 2 and m s = (1400 ± A6 stat ± 53 syst ) MeV/c 2 in 
agreement with the K*(892) and K*(1430) masses, respectively. In Fig. 7 is 
shown a comparison between our results and those of [9] and [11] for this 
parametrization. 

Using the recently proposed parametrization based on a dispersive approach, 
we obtain for the slope of the vector form factor: A + = (23.3 ± 0.5 sta t ± 
0-8 S2/s j) x 10~ 3 and for the logarithm of the scalar form factor at the Callan- 
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Treiman point: InC = (143.8 ± 8.0 sta t ± 11.2 syst ) x 10~ 3 . According to the 
model proposed in [4] the value of InC can be used to test the existence 
of RHCs by comparing it with the Standard Model predictions. Taking the 
value of \Fk+V us / F n+ Vud\ from [20] and those of / + (0)|\4s| an d \V u a\ from [21] 
we obtain for a combination of the RHCs couplings and the Callan-Treiman 
discrepancy (Act) the value: 2(es — €ns) + ^ct = —0.071 ± 0.014jva48 ± 
0.002^ eo ± 0.005 ex t, where the first error is the combination in quadrature 
of the statistical and systematical uncertainties, the second one refers to the 
uncertainties related to the approximations used to replace the dispersion 
integrals and the last one is due to the external experimental input. 
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Fig. 6. 1 a contour plots in the plane A+-A+ showing the NA48 results together with 
those of [9,7,11] for the quadratic fits of the andK c ^ decays. The dot represents 
A + and A + as obtained with a Taylor expansion of the pole parametrization. The 
ISTRA+ results have been multiplied by the ratio (m n +/m n o) 2 . 
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Fig. 7. NA48 pole fits results, together with those of [9] and [11], /orK^ and K e 3 de- 
cays. K e 3 results appear as vertical bands on the my^ms plane; the ellipses shown 
for K M 3 results are the 68% C.L. contour plots. 
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